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The neutral B-meson pair produced at the T(4S') should exhibit a non-local correlation of the type 
discussed by Einstein, Podolski, and Rosen. We measure this correlation using the time-dependent 
flavour asymmetry of semileptonic decays, which we compare with predictions from quantum 
mechanics and two local realistic models. The data are consistent with quantum mechanics, and 
inconsistent with the other models. Assuming that some B pairs disentangle to produce B" and 
with definite flavour, we find a decoherent fraction of 0.029 ±0.057, consistent with no decoherence. 

PACS numbers: 03.65.Ud, 03.65.Yz, 13.25.Hw 



The concept of entangled states, which cannot be de- 
scribed as product states of their parts, was born with 
Quantum Mechanics (QM). In 1935 Einstein, Podolski 
and Rosen (EPR) considered such a pair of particles and 
concluded that QM cannot be a "complete" theory 
this suggests that additional ( "hidden" ) variables are re- 
quired. In 1964 J. S. Bell showed that QM can violate a 
certain inequality, which is (by contrast) satisfied by all 
local hidden variable models [21 ■ Many experiments have 
since been performed and found excellent agreement with 
the prediction of QM (although no "loophole" -free Bell 
test has yet been performed) [3|. Most of these studies 
have used pairs of optical photons; it is also interesting 
to test EPR correlations in massive systems Q at much 
higher energies 0]. In this Letter, we present a study 
of EPR correlation in the flavour of iJ-meson pairs pro- 
duced at the T(45'). Contrary to the analysis presented 
in 0, and as discussed in the literature a Bell in- 
equality test cannot be performed in this system due to 
the rapid decrease in time of the i3-meson amplitudes, 
and the passive character of the flavour measurement, 
via reconstruction of i3-meson decay products. In- 
stead, we compare the data with predictions from QM 
and other models. Related studies have been performed 
in the iiT-meson system [1, Q to test decoherence (loj 
effects; T{AS) B^B^ data have also been analyzed, 
but using time- integrated information only [ll| . Here, we 
use information on reconstructed _B-meson decay times to 
test both decoherence and the Pompili-Selleri model [l^l , 
which re pres ents a range of possible local hidden- variable 
theories jl3j . 

The wavefunction of a B^B^ pair from T(4S') decay 
is analogous to that of photons in a spin-singlet state 
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(1) 



Decays occurring at the same proper time are fully cor- 
related: the flavour-specific decay of one meson fixes the 
(previously undetermined) flavour B^ / B^ of the other 
meson. Given ([1]), the time-dependent rate for decay 
into two flavour-specific states Ri = e~^*/'^fl°/(4rso){l± 
cos(A77idAt)} for opposite flavour {B^W; +, i = OF) 
and same flavour or WW; -, i = SF) decays. 

At = 1^1 — ^2! is the proper-time difference of the decays, 
and Arrid the mass difference between the two B^-B^ 
mass eigenstates. We have assumed a lifetime difference 
ATd = and neglected the effects of CP violation in 
mixing, which are 0(10"'') or less. 

Thus in QM the time-dependent asymmetry A{A.t) = 
{RoF — Rsf) / {RoF + Rsf) = cos(AmrfAt), is a func- 
tion of At but not the individual times ti.2. This is 
a manifestation of entanglement. By contrast, we can 
consider Spontaneous Disentanglement (SD), an extreme 
case of decoherence, in which the i3-meson pair imme- 
diately separates into a B'^ and 5° with well-defined 
flavour, which then evolve independently [lij . The asym- 
metry becomes 

AsB{ti,t2) = cos(Amdii)cos(Amrfi2) (2) 

= i[cos(A?n<j(ti + 12)) + cos(ATO<iAi)], 

depending on ti -\- 12 in addition to At. Because of the 
large uncertainty on the T(4S') decay point, it is difficult 
to measure individual decay times ti^2'- only At is mea- 
sured in this analysis. If we first integrate the OF and 
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FIG. 1: Time-dependent asymmetry predicted by (QM) quan- 
tum mechanics and (SD) s pon taneous and immediate disen- 
tanglement of the B-pair [13, [Til, and (PS min 

to PS max I 

the range of asymmetries allowed by the Pompili-Selleri 
model Amd — 0.507 ps^^ is assumed [23 |. 



SF distributions keeping At constant we obtain the asym- 
metry curve shown in Fig. [1] which differs significantly 
from the simple cosine term due to QM (also shown). 

In the model of Pompili and Selleri (PS) [l^, each B 
has well-defined flavour, or B^ , and mass, correspond- 
ing to the heavy and light B^-B'^ eigenstates. There are 
thus four basic states: B^, B^, B^, B^. At equal times 
At = 0, the S-mesons in a pair have opposite values of 
both mass and flavour; mass values are stable, but the 
flavour can change, simultaneously for the two mesons. 
There are no other assumptions, except a requirement 
that QM predictions for uncorrelated S-decays are re- 
produced. This rather general scheme includes a range of 
possible local realistic models, and allows time-dependent 
asymmetries to lie within the bounds 

A5?s^'^(ti,<2) = 1 - |{1 - cos(AmdAt)} cos(Amdtmi„) 

+ sin(AmdAt) sin(ATOdtmin)|, and (3) 
Ap^"(ti, <2) = 1 - min(2 + 2 - f), where (4) 

'5 = {1 + COs(ATOc( At)} COs(ATOc;tmi,i) 

— sin(Am^At) sin(Amrftn-,in). (5) 

Note the additional = min(ti,t2) dependence. Af- 
ter integration for fixed values of At we obtain the asym- 
metry curves PSmax and PSmin shown in Fig. [1] 

To determine the asymmetry, we use 152 x 10^ BB 
pairs collected by the Belle detector at the T(4S') res- 
onance at the KEKB as ymm etric-energy (3.5 GeV on 



8.0 GeV) e+e- collider [l7l. The Belle detector [1 



is a large-solid-anglc spectrometer consisting of a sili- 
con vertex detector (SVD), central drift chamber (CDC), 
aerogel Cherenkov counters (ACC), time-of- flight coun- 
ters (TOF), and a CsI(Tl) electromagnetic calorimeter 
(ECL) inside a 1.5T superconducting solenoid. The flux 
return is instrumented to detect and identify muons 
(KLM). The T(4S') is produced with p-f = 0.425 close to 
the z axis (defined as anti-parallel to the positron beam 



line). As the B momentum is low in the T(4S') center- 
of-mass system (CMS), At can be determined from the 
z-displacement of B-decay vertices: At w Az/p-fc. 

We use an event selection similar to that of a previ- 
ous Belle analysis [3, [IS] , but optimised for theoretical 
model discrimination; in particular we use more strin- 
gent criteria on the flavour tag purity than the previ- 
ous analysis. To enable direct comparison of the result 
with different models, we subtract both background and 
mistagged-flavour events from the data, and then correct 
for detector effects by deconvolution. 

We determine the flavour of one neutral B by recon- 
structing the decay B° -> D*~(.+v, with D*' ^ D"Tr~ 
and 1?° ii'+TT" (tt^) or /•C+Tr^Tr+TT" (charge-conjugate 
modes are included throughout this Letter). Charged 
particles (except the "slow pion" tTs) are chosen from 
tracks with associated SVD hits and radial impact pa- 
rameter dr < 0.2 cm, and required to satisfy kaon/pion 
identification criteria usin g co mbined TOF, ACC and 
CDC [dE/dx) information [2l|. 7r° 77 candidates are 



selected with \M. 
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Troj < llMeV/c^ and momenta 
p-^a > 0.2GeV/c; the photons must have energies > 
80MeV. We select D° candidates with (Mif„^ -to/jo) e 
[-13,13]MeV/c2 for KTT{'Kn) and [-37, 23] MeV/c^ for 
KuTT^ . A D* candidate is formed by constraining a 
and slow pion (having opposite charge to the lep- 
ton) to a common vertex. We require a mass differ- 
ence A'/diff = Mku^^, - Mku^ e [144.4, 146.4] MeV/c^, 
and CMS momentum p*jj, < 2.6 GeV/c, consistent with 
-B-decay. Electron identification uses momentum and 
dE/dx information, ACC response, and energy deposi- 
tion in the ECL. Muon identification is based on pene- 
tration depth and matching of hits in the KLM to the 
extrapolated track. The efficiency is about 92% (84%) 
for electrons (muons) in the relevant momentum region, 
from 1.4 to 2.4 GcV/c in the CMS; hadrons pass this se- 
lection with an efficiency of 0.2% (1.1%). We require that 
the CMS angle between the D* and lepton be greater 
than 90°. From the relation Ml = {E*g - E*^,^Y - 
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\cofi{9B,D'i), where Ob,d* 



is the angle between pg and Pq,^, we can reconstruct 
coa{9B.D*e) by assuming a vanishing neutrino mass. We 
require \ cos{6B,D*i)\ < 1.1. The neutral B decay po- 
sition is determined by fitting the lepton track and D'^ 
trajectory to a vertex, constrained to lie in the e+e^ in- 
teraction region (smeared in the r — (f> plane to account 
for the B flight length); we require /udof < 75. 

The remaining tracks are used to determine the sec- 
ond B decay vertex and its flavour, using the method of 
Ref. 21I, [23]. Events are classified into six subsets ac- 
cording to the purity of the tag. In this analysis we use 
only leptonic tags from the highest purity subset. 

In total, 8565 events are selected (6718 OF, 1847 SF). 
A GEANT-based Monte Carlo (MC) sample assuming 
QM correlation, with five times the number of events, was 
analysed with identical criteria; its Az and D* mass dis- 
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tributions were tuned to those of the data. This sample 
was used for consistency checks, background estimates 
and subtraction, and to build deconvolution matrices. 

To compensate for the rapid fall in event rate with 
A^, the time-dependent distributions are histogrammed 
in 11 variable-size bins (Table [T|. Background subtrac- 
tion is then performed bin-by-bin; systematic errors are 
likewise determined by estimating variations in the OF 
and SF distributions, and calculating the effect on the 
asymmetry. Terms due to event selection are estimated 
by comparing data and MC distributions for each quan- 
tity, and converting discrepancies into yield variations: 
effects due to each selection are added in quadrature. 
Estimation of the remaining terms is described below. 

Four types of background events have been considered: 
e^e^ — > qq continuum, non-D* events, wrong Z)*-lepton 
combinations, and D**'^iv events. Off-resonance 

data (8.3 fb~^) were used to estimate the continuum 
background, which was found to be negligible. 

The background to the D° sample, and misassigned 
slow pions, produce a background under the D* peak in 
Afdiff- As a correction, we subtract 126 ±6 (54 ±4) such 
OF (SF) events based on scaled yields from the sideband 
Mdiff e [156.0, 164.0] MeV/c^. The corresponding sys- 
tematic uncertainty is estimated by considering statisti- 
cal fluctuations, and moving cuts by ±0.02 MeV/c^ (the 
estimated mis-calibration in Mdiff). Alternate sidebands 
[152.0, 156.0] and [164.0, 168.0] McV/c^ are also used: the 
difference from default results (consistent with statistical 
fluctuations) is conservatively included in the systematic 
error. 

The wrong £)*-lepton combination background is 
mainly due to the combination of a D* from one B with 
a true lepton from the other B, with a smaller fraction 
due to misidentified leptons, and from charm decay. To 
estimate this background, for each selected lepton which 
forms a CMS angle to the D* less than 90°, we reverse 
its CMS momentum labelling the modified lepton , and 
require \cos{9b.d*i')\ < l-l- This procedure, intended 
to reject correlated D*£ pairs while selecting events with 
no angular correlation, has been validated on MC events 
where true B^ D*^t^Xv combinations have been ex- 
cluded. (The correlated background from charm decays 
is negligible.) We obtain 78 ± 9 OF and 237 ± 15 SF 
events, which are then subtracted. Contributions to the 
systematic error are obtained by considering the statisti- 
cal fluctuations and by moving cuts by ±0.1 to account 
for possible data-MC discrepancies. 

After these subtractions, three main types of events 
remain: B° D*-£+v, the signal; B° D**-i+v, 
which we retain because it undergoes mixing; and B^ —t 
D**^t^v background. MC shapes for the signal and the 
sum of the D** channels are used in a two-parameter 
flt to the cos{9B.D*t) distribution to find the total D** 
contribution {'X^ /udof = 56/46), and its B^ component 
is then estimated using MC fractions. We find 255.5 ± 



16.0 events (254.0 OF and 1.5 SF), which we subtract 
from the data. The systematic uncertainty is estimated 
by adding in quadrature the flt error (6%) and variations 
obtained by moving the flt region (3%) and changing to a 
single parameter flt with forced normalisation (2%). We 
also assign a 20% uncertainty on the ratio of branching 
fractions of B^ D**-£+u to B+ D**°£+v. 

We correct for wrong flavour assignments using OF 
and SF distributions from wrongly-tagged MC events. 
The mistag fraction 0.015±0.001 (stat) is consistent with 
that in data [13]; we assign a systematic error of ±0.005. 

Remaining reconstruction effects (e.g. resolution in At, 
selection efficiency) are corrected by deconvolution, treat- 
ing the SF and OF distributions separately. The method 
is based on Deconvolution with Singular Value Decom- 
position (DSVD) [2^; 11x11 response matrices are built 
separately for SF and OF events, using MC D*£i' events 
indexed by generated and reconstructed At values. The 
procedure has been optimised by a Toy Monte Carlo 
(TMC) technique where sets of several hundred simulated 
experiments are generated with data and MC samples 
identical in size to those of the real experiment, but as- 
suming different true asymmetries Aqm, ^sd, and Apg'^. 
In particular the following points have been studied: 

(1) The effective matrix rank was reduced from 11 to 5 
(6) for the OF (SF) sample, to minimize the total error. 
(The statistical precision of some singular values is poor.) 

(2) The MC events used to fill the response matrix, 
and provide an a priori to the regularization algorithm, 
introduce a potential bias: e.g. the first At bin contains 
few SF events for QM, but is well-populated for SD. We 
therefore replace SF and OF samples with mixtures SF + 
o X OF and OF + s x SF, choosing s = o = 0.2 to minimize 
systematic effects; the exact values are not critical. 

(3) After DSVD, measured differences from input val- 
ues are averaged (over QM, SD, and PS) and subtracted 
bin-by-bin from the asymmetry, to reduce the potential 
bias against any one model. The maximal absolute devi- 
ation of the corrected distribution from the three models 
is assigned as the systematic error in each At bin. 

(4) A 46 fXTa Gaussian smearing term, inferred from the 
difference between MC and data vertex-fit errors, is used 
to tune the MC Az distribution to the data. (The aver- 
age Az resolution is « 100 /^m). This term was varied by 
its ±35 /xm uncertainty, and the resulting bin-by-bin dif- 
ference in the asymmetry taken as the systematic error. 

Terms from (3) and (4) are added in quadrature to 
give the total systematic error due to deconvolution. We 
test the consistency of the method by fitting the B^ decay 
time distribution (summing OF and SF samples), leaving 
the B° lifetime as a free parameter. We obtain 1.532 ± 
0.017(stat) ps, consistent with the world average 24 1. 
We also repeat the deconvolution procedure using events 
with better vertex fit quality, and hence more precise At 
values: consistent results are obtained. 

The final results, which may be directly compared with 
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TABLE I: Time-dependent asymmetry in At bins, corrected for experimental effects, witli statistical and systematic uncertain- 
ties. Contributions from event selection, background subtraction, wrong tag correction, and deconvolution are also shown. 



At bin 


window [ps] 


A and total error 


statistical error 


total 


event sel. 


Systematic 
bkgd sub. 


errors 
wrong tags 


deconvolution 


1 


0.0 - 


0.5 


1.013 ± 0.028 


0.020 


0.019 


0.005 


0.006 


0.010 


0.014 


2 


0.5 - 


1.0 


0.916 ± 0.022 


0.015 


0.016 


0.006 


0.007 


0.010 


0.009 


3 


1.0 - 


2.0 


0.699 ± 0.038 


0.029 


0.024 


0.013 


0.005 


0.009 


0.017 


4 


2.0 - 


3.0 


0.339 ± 0.056 


0.047 


0.031 


0.008 


0.005 


0.007 


0.029 


5 


3.0 - 


4.0 


-0.136 ±0.075 


0.060 


0.045 


0.009 


0.009 


0.007 


0.042 


6 


4.0 - 


5.0 


-0.634 ± 0.084 


0.062 


0.057 


0.021 


0.014 


0.013 


0.049 


7 


5.0 - 


6.0 


-0.961 ±0.077 


0.060 


0.048 


0.020 


0.017 


0.012 


0.038 


8 


6.0 - 


7.0 


-0.974 ± 0.080 


0.060 


0.053 


0.034 


0.025 


0.020 


0.025 


9 


7.0 - 


9.0 


-0.675 ±0.109 


0.092 


0.058 


0.041 


0.027 


0.022 


0.022 


10 


9.0 - 


13.0 


0.089 ±0.193 


0.161 


0.107 


0.067 


0.063 


0.038 


0.039 


11 


13.0 - 


20.0 


0.243 ± 0.435 


0.240 


0.363 


0.145 


0.226 


0.080 


0.231 




At [ps] At [ps] At [ps] 

FIG. 2: Bottom: time-dependent flavour asymmetry (crosses) and the results of weighted least-squares fits to the (left to right) 
QM, SD, and PS models (rectangles, showing ±1(7 errors on ArUd)- Top: differences A = Adata — A^-nodci in each bin, divided 
by the total experimental error atot. Bins where Apg" < Adata < Apf have been assigned a null deviation: see the text. 



theoretical models, are shown in Table HI addition in 
quadrature is used to combine the various error terms. 

We perform weighted least-squares fits to A{At), in- 
cluding a term taking the world-average Am^ into ac- 
count. To avoid bias we discard BaBar and Belle mea- 
surements, which assume QM correlations: this yields 
(Amd) = (0.496 ± 0.014) ps-i jH]. 

In fits to the QM, SD, and PS predictions, we obtain 
Amd = 0.501 ± 0.009, 0.419 ± 0.008, and 0.447 ± 0.010 
ps""'^ with of 5.2, 174, and 31.3 respectively, for eleven 
degrees of freedom: see Fig. [2] The data favour QM over 
the SD model at 13cr, and QM over the PS model at 
5.1(7 [26| . As noted above, CP violation in mixing can 
be neglected. Introducing a lifetime difference = 
0.009 ± 0.037 [iil has a neghgible effect on the fit. As a 
consistency check, the time-dependent asymmetry before 
deconvolution is compared to MC predictions for QM and 
{via reweighting) the SD and PS models: QM is strongly 
favoured. 

Following other phenomcnological studies of decohcr- 
ence {e.g. Ref. 0) we also fit the data with the function 



(1 - Cbo5o)Aqm + Cboro^SB- this is e_quivalent to mod- 
ifying the interference term in the B'^-B^ basis, or to as- 
suming that only a fraction of the neutral B pairs from 
T(4S') decays disentangle immediately into a and a 
W. We find QgOBO = 0.029 ± 0.057, consistent with no 
decoherence. 

In summary, we have analysed neutral B pairs pro- 
duced by T(4S') decay, determined the time-dependent 
asymmetry due to flavour oscillations, and corrected for 
experimental effects by deconvolution: the results can be 
directly compared to theoretical models. Any local re- 
alistic model including the assumptions of Pompili and 
Selleri is strongly disfavoured compared to quantum me- 
chanics. Immediate disentanglement, in which definite- 
fiavour _B" and B^ evolve independently, is ruled out; if a 
fraction of B-pairs is assumed to decay incoherently, we 
find a decoherent fraction consistent with zero. 
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